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1. Introduction 

The mean field a^— dynamo of magnetoliydrodynamics (MHD) ^ |2 El plays a 
similarly paradigmatic role in MHD dynamo theory like the harmonic oscillator in 
quantum mechanics. In its kinematic regime this dynamo is described by a linear 
induction equation for the magnetic field. For spherically symmetric a— profiles 
a(r) the vector of the magnetic field can be decomposed into poloidal and toroidal 
components and expanded in spherical harmonics. After additional time separation, 
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the induction equation reduces to a set of ^— decoupled boundary eigenvalue problems 

2laU = Au, u(r \ 0) = u(l) = (1) 
for matrix differential operators 



with HE] 



Ai,c. ■■= -dra{r)dr + a(r)^^^-til = a{r)Ai - a'(r)9, . (3) 

The boundary conditions in are idealized ones and formally coincide with those 
for dynamos in a high conductivity limit of the dynamo maintaining fluid/plasma 
We will restrict our subsequent considerations to this case and assume a domain 

I?(2t„) = {u e W = L2(0, 1) ® ^2(0, 1)1 u(r \ 0) = u(l) = o} (4) 

in the Hilbert space {Ti, (., .)). The a— profile a(r) is a smooth real function C^(0, 1) 3 
a{r) : (0,1) M and plays the role of the potential in dynamo models. 
Due to the fundamental symmetry of its differential expression ^ El j 

21, = mU, = ( ? i ) ' (5) 

the operator 21q is a symmetric operator in a Krein space (/C, [.,.]) [71 151 lUl [TUl [TT] 

with indefinite inner product [.,.] = (J-t) E^nd for the chosen domain ^ it is also 
selfadjoint in this space 

[^aX,y] = [x^^ay], x,y e K.. (6) 

Below we analyze the spectrum of the operator 21^ in the vicinity of constant 
a— profiles — analytically with the help of a perturbation theory as well as numerically 
with a Galerkin approximation. We obtain a pronounced a— resonance pattern 
in the occurring deformations of the spectral mesh as well as in the unfolding of 
the semi-simple (diabolical J21) degeneration points which form the nodes of this 
mesh. Additionally, we develop a Fourier component based estimation technique for 
(/ — 0)— models which allows to obtain the critical a— profiles at which the eigenvalues 
enter the right spectral half-plane with non- vanishing imaginary components (at which 
overcritical oscillatory dynamo regimes form). 

2. Basis properties of the eigenfunctions in case of constant a profiles 

For constant a— profiles a{r) = =const, r E [0, 1), the operator matrix (j^J takes 
the simple form 

so that the two-component eigenfunctions u(r) can be easily derived with the help of 
an ansatz 

u„(r) = ^0 ®^i«(r) eC2®L2(0,l), (8) 
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Figure 1. Eigenvalues A„ (oq ^ 0) for 



(a) and i = 1 (b). The 



spectral branches intersect at semi-simple degeneration points (diabolical points) 
of algebraic and geometric multiplicity two. 



where Ci , C2 G C are constants to be determined and u„ are the eigenf unctions of the 
operator Ai 

Aun^PnUn, u(r \ 0) = u(r = 1) = . (9) 

These eigenfunctions m„ are Riccati-Bessel functions ^21 

V2 



(10) 



and we ortho-normahzcd them as^) 

{Ura,Un) ^ 5„in, |lw„|l=l. (11) 

Accordingly, the spectrum of the operator Ai consists of simple positive definite 
eigenvalues p„ > — the squares of Bessel function zeros 

J,+ i(Vp^)=0, Q<^i<^<---. (12) 

Spectrum and eigenvectors of the operator matrix SIqq follow from 0, (jHJ as 

(13) 

(14) 



and 



A„ = A„ (ao) = -pn ± aoVP" n^I/ 
^2 (g) 2,2(0,1), 



1 



and correspond to Krein space states of positive and negative type 



,± ,.±1 



,u^] = ±2^5^n, [u±,uj]=0, u±e/C±c/C. (15) 

The branches of the spectrum are real- valued linear functions of the parameter 
ao with slopes rL^fp^ and form a mesh-like structure in the (ao, 5RA)— plane, as depicted 
in Fig. In order to calculate the intersection points of the spectral branches (the 
nodes of the spectral mesh), we introduce the following convenient notation 

A^j = - Pn + eao ■, £ = ± 
1 

^\fP'n 

^) See [Appendix A| for Riccati-Bessel functions and related orthogonality conditions. 



U,, = 



(16) 
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which aUows us to treat positive and negative Krein space states in a unified way. 
Two branches A^, Afj with n ^ m intersect at a point [a^, Aq) when 

\£ _ \i5 



"0 



Pn-Pr, 



otQ ■■= OLQ = e^/^+5^/p^ 



(17) 



and hence 



Eqs. (|17|l and l|18() imply that spectral branches of different type 5 ^ e intersect 
for both signs of ap at Aq < 0. In contrast, intersections at Ag > are induced by 
spectral branches of positive type when ao > 0, and by spectral branches of negative 
type when ao < 0. 

According to equation H14|) the double eigenvalue Aq possesses the two distinct 
eigenvectors and u^: 



Consequently, the intersection points given by (|17|l correspond to double eigenvalues 
p8|l with two linearly independent eigenvectors (|19|l . i.e. they are semi-simple eigen- 
values or diabolical points |12[ll4l[THlll6| of algebraic and geometric multiplicity two. 



3. Unfolding diabolical points by perturbations of the a-profile 

Let us assume that the operator SIqjj for ao — has a semi-simple double eigenvalue 
Aq with eigenvectors and determined by equations (|18|l and (|19|l . Consider a 
perturbation of the a-profile of the form 

a{r) = < + Aa(r) = < + e^{r) . (20) 

Then, the perturbed operator is given by 

2l„=2laS+ef i^A^ u^f, "^STM + e'B (21) 

\^ ip[r)Ai — if [r)Or J ° 

and the eigenvalue problem can be expanded in terms of the small parameter e 

miTniiini as 

(2l„g + £»)« + eul + . . .) = (A;^ + eAi + . . .)« + eui + . . .) . (22) 

Here, Uq is an eigenvector of the unperturbed operator SIqjj , corresponds to the 
eigenvalue Aq and, hence, has to be a linear combination of and 

< = 7i< + 72U^ e span«, u^) C /C . (23) 
A comparison of the coefficients at the same powers of e yields up to first order in e 
St.gUQ- - K< , (24) 
2t„.ui + Q3< = A;;ui + Ai< . (25) 
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The first of these equations is satisfied identically, whereas the second one can be most 
conveniently analyzed by projecting it with the help of the Krein space inner product 
[., .] onto the two-dimensional subspace^) span(u^,u^) C JC 

- AiK,<], [S<,uf„] = AiK,ui] . (26) 

Using H23|) in (|26|l yields a closed system of defining equations for the first order 
spectral perturbation Ai and the coefficients 71 and 72 

, <] - Ai K , <] [ml , <] \ / 71 



0, (27) 



i.e. in first order approximation the perturbation e5B defines the spectral shift Ai and 
lifts the directional degeneration 123|l of the zeroth order eigenvectors Uq by fixing 
two rays in the subspace span(u^,uf„) C /C (a standard effect known also from the 
perturbation theory of degenerate quantum mechanical systems |17j'l. 

In our subsequent considerations of the system 1)27(1 we will need different 
explicit representations of the matrix elements containing 5B. Partial integration and 
substitution of the relation dfun = [—pn + ^(^ + l)/?'^]w„ give these representations as 

[*8U^, = (p [{pm + e5\/PnPrn)UmUn + + Um"n] ^r (28) 











[ed^PnPm H 5 UmUn + Uj^U, 



dr. (29) 



The symmetry properties of Eq. ((29|l and its implication 

[»ul,<] = [S<,uf;j (30) 

are a natural consequence of the Krein space self-adjointness of the perturbation 
operator 03 and the real-valuedness of the eigenvectors u^, u^. 

From H27|l and (|30|l we obtain the following defining equation for Ai 

which with the Krein space norm p5|) reduces to 

^ V '^y/P^ '^^fp^ ) ^^JPriPrn 

This quadratic equation is of the type Af — Ai(ai-|-a2)-l-aia2— e(5fe^/4 = 0, ai,a2,6e M 
and its solutions Ai_± ~ \{a\ + a-i) ± ^ (ai — a-i)^ + e(5fo^]/2 are real- valued for e — & 
and complex for e 7^ 5 n {a\ — 02)^ < 6^ (in the present first order approximation^)). 
We see that they show the typical Krein space behavior. Intersections of spectral 
branches corresponding to Krein-space states of the same type (e = (5) induce no real- 
to-complex transitions in the spectrum (they are weak interactions in the sense of il8j ). 
In contrast, intersections of spectral branches corresponding to states of different types 
(e 7^ S) may in general be accompanied by real-to-complex transitions (they are strong 
interactions in the sense of JTS). The same generic behavior is implicitly present, e.g., 
in the PT-symmetric quantum mechanical (QM) models of Refs. [Bl 1^1^ 1^1^ 

^) The terms containing ui cancel due to the self-adjointness |^ of the operator %a'^ and Eqs. 
^) Higher order corrections may lead to a further reduction of the real spectral sector. 
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(see also the discussion in The unfolding of a diabolical point in a Hermitian 

QM model under PT— symmetric perturbations was explicitly demonstrated in Ref. 

As noted above, the unfolding of the diabolical points is accompanied with a 
fixing of the directions of the zeroth order eigenvectors Uq £ span(u^,u^) C /C. Using 
Ai,± in lt?7|l one finds these directions as rays Uq j. = 71, ±u^ + 72, ±u^ defined by 

71, ± ^ [^^LK] ^ ['S^f».^fn] " K„,ui]Ai,± .^g. 

72, ± ['8u^,u^J~ K,u^]Ai,± [«uf^,u^] 

— a generic result obtained, e.g., also in pi)] . 

As part of the subsequent considerations, we will apply the general technique 
H28|l . H29(l . (|32(l and H33|l for a detailed analytical study of the unfolding of diabolical 
points in concrete dynamo setups. 



4. Local deformations of the spectral mesh 

The perturbation analysis of the previous section has been restricted to a first 
order approximation — giving trustworthy analytical results for the behavior of the 
spectrum in a very close vicinity of any single diabolical point. Here we extend this 
approximation method to parameter space regions (ao— regions) containing several 
diabolical points — allowing in this way to gain a qualitative understanding of how 
perturbations of the a— profile deform the spectral mesh over such a region. For this 
purpose we extend the projection technique of the previous section from projecting 
on two-dimensional subspaces span(u^,u^) C /C to projections on TV— dimensional 
subspaces C := span(u^i^ , . . . , u^^) C /C spanned by those eigenvectors which are 
involved in the intersections over the concrete region. The method is well known from 
computational mathematics as Galerkin method, Rayleigh-Ritz method or method of 
weighted residuals 27, 28, 29, 30, 3L. 

In order to simplify notations, we pass from double-indexed eigenvalues and 
states ufj e /Ce C /C+ © /C_ C /C with {n, e) S Z+ x Z2 ~ Z* = Z — {0} to eigenvalues'') 
A„ and normalized states 0„ depending only on the single state number^) n £ Z* : 

_ r A+ for neZ+ 9I/2 i/4„ _/ ^ /C+ for n e 1+ 

■^'^ ^ \ A|;^| for n£lr ' ^ Pn '°n-<y ^-^^ ^ ^_ fQj. nelr ^'^^> 

with obvious implication 

[Ori, Om] = e« 5nm, £n ■= sigU (n) . (35) 

Furthermore, we order the index set of the vectors of the subspace C according to the 
rule ni > n2 > ■ ■ ■ > UN with N+ vectors 0„ of positive type and iV_ = N — N+ of 
negative type. 

The approximation of the eigenvalue problem ^ consists in representing the 
eigenfunction u as linear combination 

N 

u = ^CfcO„,, CfcGC (36) 

i—k 

^) We use the notation An in spite of a possible ambiguity in the case of Ai. From the concrete 
context it will be clear whether the first order perturbation Ai is considered or the spectral branch 
A„=i. 

^) Depending on the concrete context, we will subsequently use either double-indexed or single- 
index On notations for convenience. 
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over the finite set of basis functions {Vnk}k=i ^^"^ projecting") the resulting equation 

N 

J2 Ck (2l„ - A) o„, = (37) 

1=1 

onto the subspace C — span(o„j , ■ • • , <^ I^i terms of the notation 

A[a] := (iji)5=i, := [2laO„,,0„^] , 

T] := diag(e„i, ' \ _ I In+ 



' ""^ I ~In_ 



c:^{ci,...,CNy eC'+eC'- (38) 
this leads to the simple N x A^— matrix eigenvalue problem 

(AH - A/Ar)c = 0, AH:=r7"^AH (39) 

with 

det(AH - A/at) = (40) 

as defining equation for the spectral approximation''). A few comments are in order 
here. 

First, we note that the reality of the eigenvectors t)„ and of the dynamo operator 
21q (cf. Q and (|14|l ) together with the selfadjointness of 21q in the Krein space K. 
imply that the matrix A [a] is real and symmetric, A"^[a] = A [a]. The Krein space 
related fundamental symmetry 21q = JSlJ, J (see (jSl) is reflected in the structure of 
the matrix A[q;] — T7^^A[q!] as 

AH=i7A^H»7- (41) 
The involutory matrix r) = r]^^ plays the role of a metric in the complex Pontryagin 
space^) -Kk = © C^^ 3 c. In fact, it holds for u = J2iLi Ci^m , fo = X^fLi djXinj 



u, W 



N N 

^ c^dj [t)„, , 0„J = ^ SriiCdi = cFind . (42) 

i,j = l i = l 



®) Given an exact solution Ue of the eigenvalue problem (21q — Ae)Ue = 0, the use of an approximate 
test function u = Cit)„^ leads, in general, to a non-vanishing residual (error) -R(c) := (S2la — Ae)u, 

where c := (ci, . . . , cjv)"^- The Galerkin method consists in solving the eigenvalue problem in a weak 
sense over the subspace C setting [On; , R(c)\ = for each of the vectors tJn; S £ — fixing the initially 
undefined constants c;. This implies that it yields exact solutions over the subspace C and a non- 
vanishing residual over the orthogonal complement C^^^ = KQ C The quality of the approximation 
can be naturally increased by increasing the dimension A'^ = dim£. A relatively save test for avoiding 
spurious solutions in numerical studies is to compare the output for approximations with different N. 
For details on the Galerkin method we refer to CT . 

Representing a{r) as function over a suitably chosen Af— dimensional parameter space V 9 
(Pi,...,Pm)i ct{r) = a[r;pi, . . . ,pm] the determinant approximation l4Ui would allow for easy 
studies of the unfolding behavior of the diabolical points over this parameter space. For example, 
for a— profiles a[r;pi, . . . ,pm] = ')Zi'L\Pifi{''') with the parameters pi as linear scale factors over 
a set of test functions fi(r) the determinant approximation l4Ui will lead to an algebraic equation 
F(A,pi, . . . ,Pm) = of degree deg(i^) = A'^ in the spectral parameter A and the parameters pi. With 
the help of such an algebraic equation not only the unfolding of the diabolical points can be tested 
on their sensitivity with regard to changes of the functional type fi{r), but rather the investigation 
of other (higher order) types of algebraic spectral singularities will be easily feasible. For a discussion 
(similar in spirit) on third-order branch points in "PT— symmetric matrix setups we refer to 1251 . 

*) A Krein space K, = (BK— is called a Pontryagin space 7r„ when min (dim ?C+, dim/C_) = k < 
oo [TOl. 
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Second, in the limit N± — > oo the subspace C fiUs the whole Krein space /C 
so that the approximation (|36|l of the vector u tends to the exact representation 
u = ^^=-ao^n^n ovBi the Krein space basis {Vn}^=-oo- the same limit, the 
Pontryagin space tt^ tends to the Krein space S — iS-f ® 5_ 3 c with positive and 
negative type subspaces S± as sequence spaces S± — 12{Z'^)- The determinant (|^ 
becomes a Hill type determinant. The mapping U from the eigenvalue problem ^ in 
the function space /C 9 u to its equivalent representation in the sequence space 
S 3c 

U : (2tc.-A)u = ^ (AH-A)c = (43) 

is the Krein space equivalent of the well known mapping from the quantum mechanical 
Schrodinger picture to its infinite-matrix representation in the Heisenberg picture jl7| . 
In this sense, the described Galerkin method can be understood as an approximate 
solution technique based on a 'truncated Heisenberg representation' of the eigenvalue 
problem. Obviously, the method is not restricted to a^— dynamo setups, rather in its 
present form it is applicable to any other Krein-space related setup as well, like e.g. 
models of T'T— symmetric quantum mechanics. The only ingredient needed is a set of 
exactly known basis functions of an unperturbed operator. 

In the next section, we use the described Galerkin method for a rough numerical 
analysis of the deformations of the spectral mesh in the region depicted in Fig. ^ — 
leaving analytical estimates of the residual (the approximation error) to forthcoming 
work^) . 

5. The hyper-idealized s wave [l = 0) sector and its a— resonance patterns 

In this section, we consider the hyper-idealized case of zero spherical harmonics, I = 0, 
which in analogy to quantum scattering theory can be interpreted as s— wave sector. 
Due to its too high symmetry contents, this sector does not play a role in the physics 
of spherical dynamos. There are no s— wave dynamos at all j^lE]. Instead, it can 
be understood as a disk dynamo model — in the concrete case of boundary 
conditions Q , as a disk dynamo with formal boundary conditions corresponding to a 
high-conductivity limit. Due to the strong spectral similarities of models with / = 1 
and I = (visible e.g. in Fig. the study of disk dynamo models turns out very 
instructive from a technical point of view. Due to their highly simplified structure 
they allow for a detailed analytical handling and a transparent demonstration of some 
of the essential mathematical features of the dynamo models. In our concrete context, 
they will provide some basic intuitive insight into the dynamo related specifics of the 
unfolding of diabolical points. In the (/ > 1)— sectors of the spherical models, these 
specifics will re-appear in a similar but more complicated way (see section El below). 

Subsequently, we perform an analytical study of the local unfolding of the 
diabolical points (along the lines of Section O that we supplement by numerical 
Galerkin results on the deformation of the spectral mesh. 

Let 1^0. Then the differential expression of the operator Ai-q reads simply 
Ai^Q = —d^ and Ai—q has ortho-normalized eigenfunctions u„{r), m„(0) = it„(l) = 
and eigenvalues p„ 

Un = V2sin(n7rr), p„ = (ttu)^ . (44) 

^) Collaborative work "An operator model for the MHD a^— dynamo" together with H. Langer 
and C. Tretter, (in preparation). 
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The eigenvalues of the matrix differential operator S^Iq-^ are given as 

= -(7m)^ + eaoTrn (45) 
and the corresponding eigenvectors yield Krein space inner products 

[Un,ui^] = e2nn6^sSnin (46) 
and perturbation terms 



[<Bui,<] =2^2^71 / (^(r) cos [((5' 
Jo 



m — en)T:r] dr . (47) 

According to p7|l . the diabolical points are located at points 

= TT{en + 5m) =: ttM, Aq = eSn'^mn (48) 

and form a periodic vertical line structure in the (aoi SRA)— plane ^''). 

Inspection of the defining equation (|31|l for the first order spectral perturbations 
Ai shows that this equation is invariant with regard to a re-scaling of u^, u^. Passing 
to the single-index notation H34|) yields 

A„ = —{nn)'^ + aoTin, D„(r) = —= ( ^ ] sin(7r|n|r'), n £ Z* (49) 



n 



TTU 



with [o„, Om] = Sn^nrn and the following convenient representation of the perturbation 
terms (|T7jl 

[Q5tJ„+j,D„] = 7r^|n(n-h j)| / ip{r) cos{j7Tr)dr . (50) 



This reduces the defining equation 1)32(1 for the spectral perturbation at the (n, n + j) 
node of the spectral mesh (the intersection point of the A„ and Xn+j branches of the 
spectrum) with coordinates 

ao =7r(2n+j) ttM, Ao=7r^i(n + j) (51) 

to 

Xj- \nT{2n + j) [ if{r)dr + 
Jo 

tf{r)dr ) — / ^p{r) cos{jTrr)dr 



+ 7r^ri(n -|- j) 

with solutions 
vr 

Ai,± = -{2n + j) J ip{r)dr 







(52) 



±^yj2^y (p(r)dr^ + ^n{n + j) ip{r) cos {jTTr)dr^ . (53) 

We observe that the strength of the complex valued unfolding of a diabolical point 
at a node with n{n + j) < is defined by the relation between its cos (j7rr)— filtered 
perturbation ^{r) cos(j7rr)(ir and its average perturbation ip{r)dr. 

In the (/ > 1)— sectors (for fixed I > 1) this hne structure is approached asymptotically in the 
|m|,|n| — > oo limit. It follows from substituting the / <C |n| — ► oo limit of the Bessel function 
zeros 1131 . y/P\n\ ~ \n\TT [1 + i/(2|n,|)] into the expression for the oq— coordinate I17i of the diabolical 
points: ag = £„yp|;q" -I- £myP|;;i[, m,neZ*. 
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A deeper insight into this pecuhar feature of the unfolding process can be gained 
by expanding the perturbation (p in Fourier components over the interval [0, 1] 



V(^) = "7 + X! cos(27rfcr) + bk sin(27rfcr)] (54) 



ao 

fe=i 

with coefficients given as oq = 2 ip{r)dr, Uk — 2 tp{r) cos(27rfcr)c?r, bk = 

2 fir) sin(27rfcr)dr. In this way the integral '/'(t') cos(j7rr)c?r in (|50|l reduces to 
components of the type 

cos{jTrr)dr = Sj,Q , 

1 1 
cos(27r/cr) cos{jiTr)dr = - {5j^2k + <5_j,2fc) , 



for j = ±2k 

sin(27rfcr') cos{j'Kr)dr — i (55) 



i-(-i)^ 2k 



for j ^ ±2k 



and we obtain the perturbation terms H50(l as 

1 - i-iy Abkk 



Qj '■= ao<5i,o + ^ 
fe=i 



4fc2 - j2 



(57) 



The defining quadratic equation for the spectral perturbation Ai takes now the form 

2 

Af-Ai|(2n + j>o + ^7i(n + j)(ag-Q2) =0 (58) 
and leads to the following very instructive representation for its solutions 



Ai.± - 4 



{2n + 3)ao±Jpal + An{n + j)Q] . (59) 



The structure of these solutions shows that the unfolding of the diabolical points 
is controlled by several, partially competing, effects. Apart from the above mentioned 
Krein space related feature of unfolding into real eigenvalues for states of the same 
Krein space type (n(n+j) > 0), and the possibility for unfolding into pairwise complex 
conjugate eigenvalues in case of states of opposite type {n{n + j) < 0) a competition 
occurs between oscillating perturbations {ak^o,bk) and homogeneous offset-shifts oq. 
In the case of vanishing offset-shifts (mean perturbations), ag = 0, any inhomogeneous 
perturbation with ak, bk ^ for some k > 1 leads for two branches with n{n -f j) < 
to a complex unfolding of the diabolical point. The strength of this complex directed 
unfolding becomes weaker when the homogeneous offset perturbation is switched on 
(ao 7^ 0). There exists a critical offset 

which separates the regions of real-valued and complex-valued unfoldings^'^) (in the 
present first-order approximation). For Qq < £io(c) ^ '^(^ + i) < a- complex- valued 

We note that for intersecting spectral branches it necessarily holds j ^ so that the offset term 
O'OSjfi in Qj cancels and the definition l6Ui of the critical offset Q^q^^j is justified. 




Figure 2. Odd harmonics of a(r) define the unfolding properties only of the 
colored diabolical points (equally colored DPs correspond to the same value of 
I J I), whereas even harmonics affect only the white (uncolored) DPs. 

unfolding occurs, whereas for Aq > flg^^j n n{n+j) < the diabohcal point unfolds real- 
valued. The special case of a critical (balanced) perturbation a§ = '^o{c} ^ < 
corresponds to an eigenvalue degeneration Ai_+ = Ai__ =: ^i{c) which, because of 
H33|l . (|34|l and 7i,+72^+ = li,~l2- ' coinciding zeroth order rays so that via 
appropriate normalization the corresponding vectors Oq ^, Oq _ can be made coinciding 
Dq _l_ = Dg _ =: Do(c)- This means that the original diabolical point splits into a pair of 
exceptional (branch) points at perturbation configurations oq = ±|ao(c)| with a Jordan 
chain consisting of the single (geometric) eigenvector (ray) 0^^^^ supplemented by an 
associated vector (algebraic eigenvector). This is in agreement with the unfolding 
scenario of diabolical points of general- type complex matrices described e.g. in 

Finally, the special case of j — — 2n is of interest. It corresponds to the intersection 
points located on the {ao = 0)— axis of the (ao, 3fJA)— plane, where the operator 
matrix is not only self-adjoint in the Krein space /C, but also in the Hilbert space 
7i. Due to the vanishing factor 2n + j these diabolical points unfold via perturbations 
X^^±{n/2)\n\^al-al^. 

Let us now consider the strength of the unfolding contributions induced by certain 
Fourier components. For this purpose, we note that the diabolical points at nodes 
(n,n + j) with the same absolute value of the index j are located on a parabolic 
curve ^^) 

X^o = \K'-^Y) = ^iM'~f). (61) 

Due to its special role, we will refer to j G Z as parabola index (Af S Z is the index 
of the vertical line in the (ao, SRA)— plane defined in (|48|) and l(SIJ). Furthermore, we 
see from the explicit structure of Qj (following from expression (|57|l l 

I ak [Sjqk + 5-,,2k) , 3 = ±2, ±4, . . . ^^^^ 

These {oo,j)— parametrized parabolic curves coincide exactly with the spectral curves of a 
(Z = 0)— sector model with physically realistic boundary conditions 36 . (For a discussion of physically 
realistic boundary conditions of spherically symmetric a^— dynamos we refer to ElEI-) 
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that cosine and sine components of a similar order \ai^\ \bk\ contribute differently 
at different nodes of the spectral mesh. It is remarkable that for all DPs with the 
same even \j/2\ (parabolas consisting of white points in Fig. the splitting of the 
corresponding double eigenvalues depends (modulo the pre-factors 4n(n + j), j^) only 
on the mean value oq of the perturbation ip{r) and its |j/2|-th cosine component 
ci\j/2\- For n{n + j) < these contributions are competing, whereas for a strictly 
real-valued unfolding n(n + j) > they enhance each other. Furthermore, we find 
from M := 2n + j that the even/odd mode properties of j imply the same properties 
for M: even (odd) modes affect the unfolding of diabolical points at even (odd) M 
only. This is also clearly visible from Fig. [5] 

In the more complicated case of odd parabola indices |j| the splitting of the 
diabolical points on the parabola (|61|l (in Fig. [21 they are marked as points of the 
same color) is governed by the competition between ao and the complete set of sine 
components bk of the perturbation (p{r). According to (|62|) . a dominant role is played 
by sine harmonics with « 4fc^, i.e. with /c± = (|j|±l)/2, such that a clear resonance 
and damping pattern occurs. Sine components with sa Ak^ are highly enhanced by 
a small denominator over the other sine components and the corresponding harmonics 
can be regarded as resonant ones. In contrast, sine contributions with modes away 
from the resonant k± = (|j| ± l)/2 are strongly damped by the denominator Ak'^ — 
and tend asymptotically to zero for — > oo. 

A first order approximation based insight into this asymptotical behavior can be 
gained from the positions of the exceptional points (EPs) in the (ao, 3^^^)— plane (the 
points where real-to-complex transitions occur). For this purpose, we use a ID-lattice 
type parametrization for ao in form of ag = ctg + Aao = ttM -|- Auq, M G Z and 
switch to a setting with e = 1 and Aa{r) := 0(r). Interpreting Aao as perturbation 
of a configuration with ap = ttM allows us to relate Aag to the Fourier coefficient 
ao = 2Aao and to estimate the ao— positions of the EPs Aae{M,j) relative to 
their corresponding diabolical points located on the line = ttM. Via H6Uf) and 
n = (M — j)/2 we get for given Qj 

[A«e(M,,)]^ . -^Q? • (63) 

For diabolical points in the lower (ao, 5RA)— half-plane it holds |j| > |A/|-|-2 and (|63|1 is 
well defined. In the case of cosine perturbations Ok cos(27rfcr), l|()2|l implies that only a 
single diabolical point per M — 0, ±2, ±4, . . . unfolds, whereas for sine perturbations 
bk sin(27rfcr) it leads to a countably infinite number of unfolding diabolical points per 
M = ±1,±2,.... Substituting into igSI) one obtains the KP,k^ f oo 
asymptotics of the EP positions as 

|Aa,(Af,j)|«^, (64) 

i.e. for increasing the distance of the EPs from the DPs is tending to zero. 
Conversely, 1)64(1 may be used to give an estimate for the number of complex eigenvalues 
for a given ao ~ nM + Aao close^'^) to a diabolical point line at ao ~ nM, \ Aao\ <C tt: 
j2«4|6fc|fc/(^|Aao|). 

Above, we arrived at the conclusion that both types of nodes (even ones 
and odd ones) show a similar collective behavior along the parabolic curves (|61|l 

The first-order approximation 1631 leads to a third-order polynomial F{Aae, j'^) = in j^. Its 
solutions are of limited meaning because of possible contributions from higher-order perturbation 
terms. 
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of fixed j^, responding on some specific a— perturbation harmonics in a resonant 
way. Hence, a specific a— resonance pattern is imprinted in the spectral unfolding 
picture of the diabolical points. Let us now analyze these a— resonance patterns as 




Figure 3. Resonant deformation of the spectral mesh with resonant unfolding of 
diabolical points due to perturbations by pure harmonics. Left column top down 
Aa{r) = 2.5 cos(4-7rr), Aa(r) = 2.5 cos(67rr), Aa(r) = 2.5 cos(87rr); right column 
top down Aa{r) = 2.5 sin(47rr), Ao(r) = 2.5 sin(67rr), Ao(r) = 2.5 sin(87rr). 

imprints in the deformations of the spectral mesh. We study these deformations 
with the help of a Galerkin approximation over a 24-dimensional Krein subspace 
C = span(tii2, . . . , 0-12) G JC — which is sufficient to cover the same spectral region 
as in Fig. ^ As a— profile we choose a(r) = ao + Aa(r) with pure harmonics 
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Aa(r) = fl/j cos(27rfcr) and Aa{r) = &fc sin(27rfcr), k — 2,3,4, — — 5/2 as 
perturbations^^). The expUcit structure of the corresponding Pontryagin space related 
matrix A [a] (see pyfl ) is given in [Appendix B| and yields the spectral approximations 
depicted in Fig. 

The most striking feature of the spectral deformations is their very clearly 
pronounced resonance character along parabolas with fixed index \j\ — leaving 
spectral regions away from these resonance parabolas almost unaffected. Specifically, 
we find for cosine perturbations (depicted in the left column) that the harmonics 
k = 2, 3, 4 affect only the unfolding of diabolical points located strictly on the 
associated parabolas with index j = 2k. The effect of sine perturbations with mode 
numbers k — 2, 3, 4 is shown in the right column graphics of Fig. |3| As predicted 
by (|62|l . we find a strongly pronounced unfolding of diabolical points located on the 
parabolas with |j| = 2/c ± 1, that is for |j| = 3 and |j| — 5 (upper right picture), 
\j\ =5 and |j| = 7 (middle right), and |j| = 7 and |j| = 9 (lower right picture). The 
DPs with \j\ = 2k ± m, m > 1 are less affected and the strength of the unfolding 
quickly decreases with increasing distance m to the resonant parabolas. In addition 
to the unfolding effects predicted analytically by first order perturbation theory, the 
top and middle right pictures show additional DP unfoldings on the large— ao—end of 
the \j\ = 2k parabola. The origin of these unfoldings can be attributed to higher-order 
perturbative contributions. 

The a— resonance pattern has a simple physical interpretation. Due to the fact 
that the spectral parameter implies an e^* behavior of the corresponding field mode, 
the Of— resonance pattern shows that short scale perturbations of the a— profile (Fourier 
components with higher k) coherently affect faster decaying (more negative SRA) field 
modes than large scale perturbations with smaller k (which lead to smaller negative 
3fiA). 

Up to now we used the Galerkin method for investigations of weak perturbations 
over an ao =const background. In Fig. 01 we demonstrate that the method works 
for strong perturbations as well. We observe that, increasing the strength of the pure 
harmonic perturbations from 02,62 = 5 up to 02,62 = 30, the regions with complex 
conjugate spectral contributions grow and finally intersect each other. Additionally, 
they shift into the upper (ao, 5RA)— plane leading to overcritical oscillatory dynamo 
regimes (5RA > 0, 3A 7^ 0). An estimate of critical a— profiles, for which such a 
transition to the upper (oq, 5RA)— half-plane starts to occur, can be given within a first- 
order (linear) perturbative approximation by assuming 5ftA(Q!o) ~ for the exceptional 
point closest to the (3?A = 0)-hne. Relations (EU, (ESI and = X!^ + 3?Ai = 
Ag -I- ttMoq /4 yield this condition in terms of the Fourier components Qj of such a 
critical a— profile as 

1/2 

5RA = ^ (M^ - f) + \m {1 ^) m = 0. (65) 

This relation may be used for testing concrete a— profiles on their capability to produce 
complex eigenvalues in the right spectral half-plane (3ffA > 0,3A ^ 0). 

We restrict our present consideration of {I = 0)— models to this first numerical 
output and the analytical estimate, expecting physically more relevant results from 
extending the present methods to models with physically realistic boundary conditions. 

^'*) The amplitudes 0^,6^ have been chosen as large as = bf^ = 5/2 > 1 in order to clearly 
demonstrate the unfolding pattern. 



Krein space related perturbation theory for MHD ~ dynamos 



15 




300 n 



200- Aa(r)=5sin(47tr; 

100- 




ReX 





ReX 

Aa(r)=20cos(47ir) 




300 n 

200 
100 

0" 
-100 
-200 
-300 







- Aa(r)=20sin(47rr) 




■\ 10 ~~ — 







300 
200 
100- 



Re;. 

Aa(r)=30cos(47rr; 





Figure 4. Resonant deformation of the spectral mesh under strong and ultra- 
strong perturbations of pure harmonic type: Aa(r) = 02 cos(47rr) (left column) 
and Aa(r) = 62 sin(47rr) (right column) with 02,^2 = 5,10,20,30. Strong 
perturbations lead not only to large spectral regions with complex conjugate 
eigenvalues, but also to a lifting of the corresponding complex branches into the 
upper (aoi 5RA)— plane, i.e. they lead to overcritical dynamo regimes (J?A > 0) of 
oscillatory type (3A ^ 0). 
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In the next section we present some first few results on the (/ > 1)— sectors of the 
a ^—dynamo model. 



6. Numerical techniques and examples for the {I > 1) sectors 

In this section, we reshape the general results of section |3 in a form suitable for 
numerical investigations and demonstrate them on a first concrete model from the 
(/ — 1)— sector. 

We start by representing the perturbation terms [*Bu^,u^] from Eqs. (|28|l . (|29|l 

as 

»i 



Jo 



£0\JPmPn ^ 2 ) "™"« + "m^^n 



dr 



Jo 



(66) 



or equivalently 



£m£n^JP\m\P\n\ H ) U\m\U\n\ 



dr 



(67) 



The functions are symmetric g: 
derivative^^) 



/ 

mn 



9n 



and from the Frechet (functional) 
V^[Q5o„„o„] =5^„„(r) (68) 



we find that they can be naturally interpreted as components of the perturbation 
gradient in the Krein space JC. Their explicit representation in terms of Bessel 
functions is given in [Appendix C| Representation H66|l may prove especially useful 
for the optimization of a-profiles with regard to given constraints or experimental 
requirements. 

First, we note that in terms of the gradient functions .9mn('") the defining equation 
H32|l for the first order spectral perturbations Ai reduces to 

^1 - -^1 / V {sg'nn + SgLrn) dr + 
10 



with solutions 



^9mmdr 



1 X 2' 

mn 



(69) 



f[£9\ 



S9'„ 



dr 







4^ 


Jo 





2 




2 




^9lnm) dr 


+ Ae6 


/ 'PgLndr 
Jo 




(70) 



The Frechet derivative \7f{x) of a function f{x) over an open set X 9 a; of a Banach space F D X 
is defined as f(x + u)-f{x) = V f {x)u + o{\\u\\) I^E^EU what for the functional f{ip) := [530™, D„] 
can be reshaped as /(i/p + x) ~ fiv) = C^vfif)' x) + o(llxll)- 
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Comparison with the resuhs for the hyper-ideahzed s— wave sector (disk dynamo) 
shows that apart from the generic Krein space related behavior (no complex 
eigenvalues for intersecting spectral branches of the same type, S — e, and possible 
formation of complex eigenvalues for branches of different type, e ^ 6) we find a 
generalization of the offset and oscillation contributions for e 6 type intersections. 
In rough analogy, the role of a transition preventing offset is played by the 'diagonal' 
terms 



f {gin + 



9 mm) 



dr 



(71) 



whereas the 'off-diagonal' terms 

-1 2 



(72) 



enhance a possible transitions to complex eigenvalues, i.e. a transition occurs for 



e^6 n 



in'- 



(fdr 



< 



g^inVdr 



(73) 



v2jo / \J{) 

The condition (|73|l yields an explicit classification criterion for a— profile perturbations 
with regard to their capability to induce complex eigenvalues. It may serve as an 
efficient search tool for concrete a (r)— profiles — and transforms in this way some 
general observations of Ref. 4 into a technique of direct applicability. 





Figure 5. Functions Snrnix) f°i' ' 
cases £(5 = 1 (b) and sS = —1 (c). 



1, n = m (a) and n = 1, m = 2,3 in the 



The functions glmi^) ^'^^ ' ~ 1 different values of n, m and eS are shown in 
Fig ||2Jl. One clearly sees that the 'diagonal' functions with m = n are always non- 
negative and the graphics of any two of them, (7^„(r), gmm(^) with m ^ n intersect 
only marginally (see Fig. 0). Their sums g^mi''') +ffmm(^) strictly sign-preserving 
functions so that they act as averaging integration kernels. In contrast, the 'off- 
diagonal' functions 5^m(r) with m ^ n show strong sign changes and in this way 
they act as filter kernels. Hence, the main qualitative roles of the 'diagonal' and 'off- 
diagonal' functions as 'offset' and 'oscillation filter' functions remain preserved also 
for the (/ = 1)— sector. Their subtle interplay is crucial for the a-profile to act as 
generator of complex eigenvalues. 

As in the {I — 0)— sector, spectral deformations over a certain parameter space 
region can be studied numerically. For a— profiles a{r) — a{)-\-l^a{r) the corresponding 
approximation matrix of the Galerkin method reads simply 

Amn[a] ^ KrAnn + s-ni I Aa{r)gl^^{r)dr . (74) 
Jo 
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Figure 6. Gradient function differences gn„{r) — as they enhance 

the real valued unfolding of diabolical points in the upper {ao ,'SiX)— plane (left 
figure), as well as offset gradient functions (integration kernels) of the type 
glni'^) +9n-in-ii^) for n = 2, . . . , 6 (right figure). 



300 




Figure 7. Deformations of the spectral mesh with unfolding of diabolical points 
for the [l = 1)— sector of a model with a[r) = ao + 2.5 cos(67rr). 



In Fig. 13 we illustrate the method for an a— profile a{r) = ao + 2.5cos(67rr) and 
a similar approximation subspace C — span(0i2, . . . , D-12) C /C as in the previous 
section. For the {I = 1)— mode Riccati-Bessel functions a cosine perturbation is no 
longer an exact resonance mode and the deformations of the mesh are spreading over 
a broader parabola-like region. 

Finally, we note that explicit analytical considerations of the spectrum in the 
{I > 0)— sector are obstructed by the lack of simple transformation rules between 
Riccati-Bessel functions as well as of simple expressions for integrals over triple 
products of spherical Bessel functions in case of finite integration intervals. 

7. Conclusions and discussions 



In the present work, the spectral properties of spherically symmetric MHD 
a^— dynamos with idealized boundary conditions have been studied. Using the 
fundamental symmetry of the dynamo operator matrix and the solution set of 
a model with constant a— profile, a Krein space related perturbation theory as 
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well as a Galerkin technique for numerical investigations have been developed. 
As analytical result of the first-order perturbation theory we found a strongly 
pronounced a— resonance pattern in the unfolding of diabolical points. The resonance 
behavior reflects the correspondence between the characteristic length scale of 
a— perturbations and the decay rates of the coherently induced field excitations. 
The observed correlations will strongly affect the specifics of reversal processes of 
dynamo maintained magnetic fields [l^7\ 1881 and support corresponding numerical 
simulations on more realistic dynamo setups [HHI- For the (/ = 0)— sector, a Fourier 
component based estimation technique has been developed for obtaining the critical 
a— profiles at which the eigenvalues enter the right spectral half-plane with non- 
vanishing imaginary components (at which overcritical oscillatory dynamo regimes 
form). The analytical results on the perturbative unfolding of diabolical points have 
been supplemented by numerical studies of the deformations of the dynamo operator 
spectrum. The capability of the used Galerkin approach has been demonstrated in 
extending the strength of the a— perturbations from weakly perturbed regimes up 
to ultra-strong perturbations. Extensions of the presented techniques to spherically 
symmetric a^— dynamos with realistic boundary conditions as well as to models of 
PT— symmetric quantum mechanics are straight forward. 
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Appendix A. Bessel function relations 

The solutions u„(r) of the eigenvalue problem ((HJ 

AiUn=PnUn, M„ (r \ 0) = u„ (r = 1 ) = 0, 

are so called Riccati-Bessel functions ^31 

Un{r) ^ y/p^rji{y/p:^ir) 

where ji{yf^r) denote spherical Bessel functions. They can be expressed in terms of 
Bessel functions J;_|_ i (y^p^r) and for definiteness we represent the solutions u„ as 

u„(r) = NnVrJi+i{y/p^r) (A. 3) 

with Nn a normalization coefhcient. 

The orthogonality of these solutions can be easily verified with the help of well 
known Bessel function relations |13[ 1401 HT] , For rn ^ n it holds 

(A.4) 



ium,Un) = NmN„ / rJi^i{y^r)Ji^i{y/p:^r)dr 



N„,Nn- 



Pn 



Ji+i{^/%'ir)drJi+i{^/p7ir) - Ji+i{^/p^r)drJi+i{ 
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This expression vanishes because of J;_|_i(0) — Ji^^{^/fhi) — and 
\drJi^i{y/pkr)\r\Q < oo. In a similar way it holds [OlBniHT] 



(A.5) 



so 



Un\\ =(u„,u„) = ^n/ rJi^i{y/^r)dr = -N^ drJi+ii^/p^r) 



that Nn = \/2/J;+3 (^p;;;) gives (u„,u„) = 1. 



Appendix B. Explicit structure of the approximation matrix A [a] 

In the (/ = 0)— sector, the relations (|49|l . H5U|) . (|56ll lead for the approximation 
matrix A[a] = »7~^(A„„)^ „^_^, Amn = [^a'o,n,Vn] of 123) over an a-profile 
a{r) = ao + Aa(r) to the following structure 



^mn[a] = A„i(5„i„ + e,„7rv|TOn| / Aa{r) cos[(r7i - n)TTr]dr 

JO 

= [-{nmf + aonm]S„ui+£7n-^\/\rrm\Qrn-n- (B.l) 
In the case of Aa{r) — ak cos(27rfcr) this gives: 

Amn[ao, ak] = [-{■nmf + aoTrm] 5mn + £rnl^\/\'mn\ak[5m,n+2k + ^m,n-2/c) , (B.2) 

i.e. A[a] has the eigenvalues Am of the unperturbed 21q,(, on its diagonal and possesses 
two subdiagonals with non-vanishing entries a distance ±2fc aside the main diagonal. 
For perturbations Aa(r) = 6^ sin(27rfcr) one finds 

2kbk 



+ e„,y/\^\ [l-(-l)^ 



4fc2 — (to — ny 



(B.3) 



Appendix C. Explicit expressions for the gradient functions (/^^(r) 

Explicit expressions for the gradient components gj„n('') can be derived from the 
representation I|1U|) of the eigenfunctions w^" and the defining relation As result 
one obtains in (n,£„) G Z+ x Z2 notations 

+ ;2 + J/+1/2 (f-^p^) >/; + i/2 (f-v/p^) 

r(PmPn) Jl^l/2[\/P^i)Jl^l/2[y/P^) 

_^ [ry/p^Ji-1/2 {r^/p^) - lJi+1/2 (i^Vp^)] [ry/p7iJi~i/2 {r./p:;!) - lJi+1/2 {r^/p^)] 

r {pniPnt^^ Jl-1/2 (VP™) Jl~l/2 

and gj„„(r = 1) = {p,nPnt''^- 



(C.l) 
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